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Abstract 

The Weakly Interacting Massive Particles (WIMPs) are among the main candidates 
for the relic dark matter (DM). The idea of the direct DM detection relies on elastic 
spin-dependent (SD) and spin-independent (SI) interaction of WIMPs with target 
nuclei. The importance of the SD WIMP-nucleus interaction for reliable DM de- 
tection is argued. The absolute lower bound for the detection rate can naturally 
be due to SD interaction. An experiment aimed at detecting DM with sensitivity 
higher than 10~^ event/day/kg should have a non-zero-spin target. 



1 Introduction 

The lightest supersymmetric (SUSY) particle (LSP) neutralino is assumed to be the 
Weakly Interacting Massive Particle (WIMP) and the best dark matter (DM) candi- 
date. It is believed that for heavy enough nuclei this spin-independent (SI) interaction 
of DM particles with nuclei usually makes the dominant contribution to the expected 
event rate of its detection. The reason is the strong (proportional to the squared mass of 
the target nucleus) enhancement of SI WIMP-nucleus interaction. Nevertheless there are 
at least three reasons to think that SD (or axial-vector) interaction of the DM WIMPs 
with nuclei could be very important. First, contrary to the only one constraint for SUSY 
models available from the scalar WIMP-nucleus interaction, the spin WIMP-nucleus in- 
teraction supplies us with two such constraints (see for example [1] and formulas below). 
Second, one can notice [2, 3] that even with a very sensitive DM detector (say, with a 
sensitivity of 10~^ events/day/kg) which is sensitive only to the WIMP-nucleus scalar in- 
teraction (with spin-less target nuclei) one can, in principle, miss a DM signal. To safely 
avoid such a situation one should have a spin-sensitive DM detector, i.e. a detector with 
non-zero-spin target nuclei. Finally, there is a complicated nuclear spin structure, which, 
for example, characterized by the so-called long g-tail form-factor behavior. Therefore for 
heavy target nuclei and heavy WIMP the SD efficiency to detect a DM signal is much 
higher than the SI efficiency [4]. However, simultaneous study of both spin-dependent 
and spin-independent interactions of the DM particles with nuclei significantly increases 
the chance to observe the DM signal. 



1 



X 



Earth 



Figure 1: Due to the expected annual modulation signature of the event rate (1) only the 
Sun-Earth system is a proper setup for the successful direct DM detection. 

2 Two constrainsts for SUSY due to the spin 

One believes to detect directly a DM particle x via its elastic scattering on a target nucleus 
{A, Z). The nuclear recoil energy E-^ {En ~ 10~^m^ ?s few keV) is measured by a proper 
detector (Fig. 1). The differential event rate depends on the distribution of the relic DM 
particles in the solar vicinity f{v) and the cross section of LSP-nucleus elastic scattering 
[51-112]: 

^ = r- dvnv)rA,W), Ea = «V(2Af.). (1) 

dE^ Jv^in dq^ 

Here, N = M / A is the number density of target nuclei. M and A stand for the Avogadro 
number and the atomic mass in AMU, respectively. denotes the nuclear mass, f max = 
Vesc ~ 600 km/s, Wmin = {MAEj^/2fi\f^^, the DM density = 0.3 GeV-cm-^. The 
neutralino-nucleus elastic scattering cross section for spin-non-zero (J 7^ 0) nuclei is a 
sum of the coherent (spin- independent) and axial (spin-dependent) terms [4, 13, 14, 15]: 

rf^^. 2^ ^ EI-MP Si,{q') Si,{q') 



dq^' ^ TTv\2J+l) v^{2J+l) v^{2J+l) 



It is useful to separate the zero-momentum transfer cross sections and introduce the 
normahzed-to- unity (-^sdsi(O) ~ nonzero-momentum-transfer nuclear form- factors: 



'S'sD,Sl(g^) 
'S'sD,Sl(0) 



'ksM') - ^F^- (3) 



The SD structure function <S'g^(g) contains the isoscalar 5'oo, isovector and the inter- 
ference ^oi terms: 

SioiQ) = (4iSoo{q) + ajSniq) + aoaiS'oi(g). (4) 

Here the isoscalar oq = a„ + Up and isovector ai = — a„ effective coupling constants arc 
used (see (9)). For q — the nuclear SD and SI cross sections take the forms 

^si(O) = (2J+1) =^^^si(0), (5) 

c^sdIO) = (2J + 1) ^~ — J — rp\^p/ + ■ y^) 

Here, = is the reduced x-nucleus mass and iii — ul is assumed. The 

+ Ma 

dependence on effective neutralino-quark couplings Cq and Aq in the underlying (SUSY) 
theory 

^eff = (A • XTmTsX • Ql^lbQ + Cq-xX-QQ) + - (7) 

g 

and on the spin (AjP-'^)) and the mass {fjf'^^) structure of nucleons enter into these 
formulas via the zero-momentum-transfer proton and neutron SI and SD cross sections: 



alM = ^^cl ag-(O) = 12^aJ,„; (8) 

^'" = E^^./^"^ a, = i:^,A(^), a„ = i:^,A(-). (9) 

The factors A^P'"\ which parametrize the quark spin content of the nucleon, are defined 
as 2AIj^'^^s^ = {p, s\il!q'-f^^r:,ipq\p, s)(p,„). The (Sp(„-)) is the total spin of protons (neutrons) 
averaged over all A nucleons of the nucleus {A, Z): 

(S,^(„)) ^ {A\S^^^^\A) = (A| ^s;(„)|A) (10) 



The mean velocity (v) of the relic DM particles of our Galaxy is about 300 km/s = 
10~^c. For not very heavy and Ma one can use the SD matrix element in zero 
momentum transfer limit [15, 16] 

M (X {A\apSp + anSn\A) ■ s^. (11) 



Note a coupling of the spin of x, s^, to the spin carried by the protons and the neutrons. 
The uncertainties arising from electroweak and QCD scale physics are incorporated in the 
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Table 1: Zero momentum spin structure of nuclei in different models. The measured 
magnetic moments used as input are enclosed in parentheses. Prom [17]. 



{Lj = 5i/2) 


(Sp) 


(Sn) 


H (in //at) 


ISPSM, Ellis-Flores [18, 19] 


1/2 





2.793 


OGM, Engcl-Vogcl [20] 


0.46 





(2.629)exp 


EOGM {gA/gv = 1), Engel-Vogel [20] 


0.415 


-0.047 


(2.629)exp 


EOGM {gA/gv = 1-25), Engel-Vogel [20] 


0.368 


-0.001 


(2.629)exp 


SM, Pacheco-Strottman [21] 


0.441 


-0.109 




SM, Divari et al. [22] 


0.4751 


-0.0087 


2.91 


^^Na {Lj = P3/2) 


(Sp) 


(Sn) 


ji (in hn) 


ISPSM 


1/2 





3.793 


SM, Ressell-Dean [15] 


0.2477 


0.0198 


2.2196 


OGM, Ressell-Dcan [15] 


0.1566 


0.0 


(2.218)cxp 


SM, Divari ar al. [22] 


0.2477 


0.0199 


2.22 


^■^Al {Lj = A5/2) 


(Sp) 


(Sn) 


11 (in //at) 


ISPSM, Ellis-Flores [18, 19] 


1/2 





4.793 


OGM, Engel-Vogel [20] 


0.25 





(3.642)exp 


EOGM {gA/gv = 1), Engel-Vogel [20] 


0.333 


0.043 


(3.642)exp 


EOGM {gA/gv = 1-25), Engel-Vogel [20] 


0.304 


0.072 


(3.642)exp 


SM, Engel et al. [16] 


0.3430 


0.0296 


3.584 


^^Ge {Lj = G9/2) 


(Sp) 


(Sn) 


IX (in hn) 


ISPSM, Ellis-Flores [18, 19] 





0.5 


-1.913 


OGM, Engel-Vogel [20] 





0.23 


(-0.879)exp 


IBFM, lachello et al. [23] and [14] 


-0.009 


0.469 


-1.785 


IBFM (quenched), lachello et al. [23] and [14] 


-0.005 


0.245 


(-0.879)exp 


TFFS, Nikolaev-Klapdor-Kleingrothaus, [24] 





0.34 


— 


SM (smah), Resseh et al. [14] 


0.005 


0.496 


-1.468 


SM (large), Ressell et al. [14] 


0.011 


0.468 


-1.239 


SM (large, quenched), Ressell et al. [14] 


0.009 


0.372 


(-0.879)exp 


"Hybrid" SM, Dimitrov et al. [25] 


0.030 


0.378 


-0.920 


'-'l {Lj = D,„) 


(Sp) 


(Sn) 


jjL (in iin) 


ISPSM, Elhs-Florcs [19, 26] 


1/2 





4.793 


OGM, Engel-Vogel [20] 


0.07 





(2.813)exp 


IBFM, lachello et al. [23] 


0.464 


0.010 


(2.813)exp 


IBFM (quenched), lachello ct al. [23] 


0.154 


0.003 


(2.813)exp 


TFFS, Nikolaev-Klapdor-Kleingrothaus, [24] 


0.15 





— 


SM (Bonn A), Ressell-Dean [15] 


0.309 


0.075 


2.775 {2.470}eff 


SM (Nijmegen II), Ressell-Dean [15] 


0.354 


0.064 


3.150 {2.7930}eff 


^^^Xe {Lj = L>3/2) 


(Sp) 


(Sn) 


IJi (in /Xiv) 


ISPSM, Ellis-Flores [18, 19] 





-0.3 


1.148 


OGM, Engel-Vogel [20] 


0.0 


-0.18 


(0.692)cxp 


IBFM, lachello et al. [23] 


0.000 


-0.280 


(0.692)exp 


IBFM (quenched), lachello et al. [23] 


0.000 


-0.168 


(0.692)exp 


TFFS, Nikolaev-Klapdor-Kleingrothaus, [24] 




-0.186 




SM (Bonn A), Ressell-Dean [15] 


-0.009 


-0.227 


0.980 {0.637}cfr 


SM (Nijmegen II), Ressell-Dean [15] 


-0.012 


-0.217 


0.979 {0.347}eff 


QTDA, Engel [4] 


-0.041 


-0.236 


0.70 



4 



factors Op and a„. The nuclear matrix element A4 in Eq. (11) is often related to the matrix 
element of the nuclear magnetic moment, which also consists of the matrix elements of 
the total proton and neutron spin operators: 

fi = {A\g'^Sn + glLn + g;Sp + g^LplA). (12) 

The free particle (jf-factors (gyromagnetic ratios) are (in nuclear magnetons): g^ = —3.826, 
gl^ = 0, gp = 5.586, gp = 1. The nuclear magnetic moment fi is often used as a benchmark 
for the accuracy of the calculation of Sp and S„ [14, 15]. For the most interesting isotopes 
either (S^) or (S^) dominates ((S^j-p-,) ^ {^p(n)))- S^^' f*^^ example. Table 1. 
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Figure 2: Exclusion curves for the spin-dependent WIMP-proton cross section (crgj-, as 
a function of the WIMP mass). DAMA/NaI-7a(f) contours for the WIMP-proton SD 
interaction in ^^^I are obtained on the basis of the positive signature of annual signal 
modulation [27, 28]. The scattered points are calculations of [29]. 
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Figure 3: Exclusion curves for the SD WIMP-neutron cross section {cr^j) versus WIMP 
mass) . 



From Eqs. (6) one can conclude the spin observables in DM search give us TWO 
independent constraints on a SUSY model via crgj-,(0) and (t^j){0), or, equivalently, via 
and a„. These constraints are usually presented in the form of exclusion curves obtained 
with different target nuclei (Figs. 2 and 3). There is only one similar constraint from 
spin-independent DM search experiments (Eq. (5)). This presentation is a bit obsolete 
[27, 28, 29], but it allows one to compare sensitivities of different experiments. 



3 Long-tail g-behaviour due to the spin 

As becomes larger, the finite momentum transfer limit must be considered for heavier 
mass Ma nuclei. The differential SD event rate with structure function (4) has 

now the form 



/'■''"^("'727TTw^^o(.). (13) 



dq'^ m^niA J (2 J + l)f ^ 

Comparing this formula with the observed recoil spectra for different targets (Ge, Xe, 
F, Nal, etc) one can directly and simultaneously restrict both isoscalar and isovector 
neutralino-nucleon effective couplings ao,i. These constraints will impose most model- 
independent restrictions on the MSSM parameter space. Another attractive feature of 
the SD WIMP-nucleus interaction is the g-dependence of SD structure function (4). The 
ratio of SD to SI rate in the ^^Ge detector grows with the WIMP mass [2, 3]. The growth is 
much greater for heavy target isotopes like xenon. The reason is the different behavior of 
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the spin and scalar structure functions with increasing momentum transfer. For example, 
the xenon SI structure function vanishes for ^ 0.02 GeV, but the SD structure function 
is a non-zero constant in the region (Fig. 4). As noted by Engel in [4], the relatively long 
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Figure 4: The ^^^Xe structure function for a pure bino neutralino. The single-particle 
structure function is normalized to 5(0) = 1. From [15]. 

tail of the SD structure function is caused by nucleons near the Fermi surface, which do 
the bulk of the scattering. The core nucleons, which dominate the SI nuclear coupling, 
contribute much less at large q. Therefore the SD efficiency for detection of a DM signal 
is higher than the SI efficiency, especially for very heavy neutralinos. 

4 One does not miss a DM signal due to the spin 

To estimate the DM detection rate we traditionally use the so-called effective scheme of 
MSSM (effMSSM) whose parameters are defined directly at the electroweak scale, relaxing 
completely constraints following from any unification assumption (see, for example [30]- 
[36]). Our MSSM parameter space is determined by the entries of the mass matrices of 
neutralinos, charginos, Higgs bosons, sleptons and squarks. The relevant definitions can 
be found in [31]. We have included the current experimental upper limits on sparticle and 
Higgs masses from the Particle Data Group. Also, the hmits on the rare 6 — > 57 decay 
have been imposed. For each point in the MSSM parameter space (MSSM model) we 
have evaluated the relic density of the light neutralinos ^^^^o with our code [37, 38, 39] 
based on [40], taking into account all coannihilation channels with two-body final states 
that can occur between neutralinos, charginos, sleptons, stops and sbottoms. We assume 
0.1 < Q-Ji^ < 0.3 for the cosmologically interesting region and we also consider the 
WMAP reduction of the region to 0.094 < fi^/?^ < 0.129 [41, 42]. 

From Fig. 5 one sees that the SD contribution obviously dominates in the domain of 
large expected rates in the non-zero-spin germanium detector {R > 0.1 event /day/kg). 
But as soon as the total rate drops down to R < 0.01 event/day/kg or, cquivalently, 
the SI (scalar) neutrahno-proton cross section becomes smaller than 10~^-10""^'^ pb, the 
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Figure 5: Ratio of the SD (spin) event rate to the SI (scalar) event rate in the ^'^Ge isotope 
(spin= 9/2) as a function of the total (SD+SI) event rate (left) and the scalar cross section of 
the neutralino-proton interaction (right). The solid vertical lines give the expected sensitivity 
of one of the best future projects GENIUS [43] . In the region above the horizontal line the spin 
contribution dominates. 

SD interaction may produce a rather non-negligible contribution to the total event rate. 
Moreover, if the scalar cross section further decreases {a < 10^^^ pb), it becomes obvious 
that the spin contribution alone saturates the total rate and protects it from decreasing 
below R ^ 10~^-10~'^event/day/kg [44]. With only a spinless detector one can miss a 
signal caused by SD interaction. An experiment aimed at detecting dark matter with 
sensitivity higher than 10~^event/day/kg should have a non-zero-spin target. Indeed, 
while the scalar cross sections governed mostly by Higgs exchange can be rather small, the 
spin cross section cannot be arbitrarily small because the mass of the Z boson [29] , which 
makes the dominant contribution, is well defined, provided one ignores any possible fine- 
tuning cancellations. Therefore, if an experiment with sensitivity 10~^-10~^ event/day/kg 
fails to detect a dark matter signal, an experiment with higher sensitivity should have a 
non-zero-spin target and will be able to detect dark matter particles only due to the spin 
neutralino-quark interaction. 



5 Conclusion 

There are at least three reasons to think that spin-dependent interaction of the DM 
WIMPs with nuclei could be very important. First, contrary to the only one constraint 
for SUSY models available from the spin-independent WIMP-nucleus interaction, the SD 
WIMP-nucleus interaction supplies us with two such constraints. Second, for heavy target 
nuclei and heavy WIMP masses the SD efficiency to detect a DM signal is much higher 
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than the SI efficiency. Finally, the absolute lower bound for the DM detection rate can 
naturally be due to SD interaction. An experiment aimed at detecting DM with sensitivity 
higher than 10^^ event/day/kg should have a non-zero-spin target. 

The authors have a pleasure to thank Prof. H.V.Klapdor-Kleingrothaus for fruitful 
and permanent cooperation as well as the organizers of the Baldin Seminar for their 
invitation to give this pedagogical review talk whose subject is beyond the traditional 
scope of the Seminar. 
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